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Besides the traditional circuit-based model of quantum computation, several quantum algorithms 
based on a continuous-time Hamiltonian evolution have recently been introduced, including for in- 
stance continuous-time quantum walk algorithms as well as adiabatic quantum algorithms. Unfor- 
tunately, very little is known today on the behavior of these Hamiltonian algorithms in the presence 
of noise. Here, we perform a fully analytical study of the resistance to noise of these algorithms using 
perturbation theory combined with a theoretical noise model based on random matrices drawn from 
the Gaussian Orthogonal Ensemble, whose elements vary in time and form a stationary random 
process. 



I. INTRODUCTION 

There has been a growing interest recently for the con- 
cept of Hamiltonian-based quantum algorithms, as op- 
posed to the standard circuit-based paradigm of quantum 
computing. The Hamiltonian algorithms rely on the con- 
tinuous time evolution of a quantum register according to 
the Schrodinger equation, and include in particular the 
quantum search algorithms by adiabatic evolution yij or 
by continuous quantum walks 0, . While these algo- 
rithms may be translated into circuit-based algorithms so 
that they could be implemented on a "standard" quan- 
tum computer 0, Q , another possibility is to consider a 
"continuous" quantum computer specifically designed to 
run this type of algorithms [6| . For a realistic implemen- 
tation, it seems therefore crucial to investigate how well 
such a quantum computer would behave in the presence 
of noise. Until now, this question has only been addressed 
for some specific algorithms subject to some very partic- 
ular noise. For instance, Childs et al. have considered an 
adiabatic quantum algorithm for solving combinatorial 
problems p| affected by an error modeled by an extra 
term which is random but deterministically evolves in 
time [8|- While this study was purely numerical, later 
on Shenvi et al. [j| analytically analyzed the effect of a 
Markovian stochastic variable perturbing the amplitude 
of the oracle Hamiltonian in the specific case of the analog 
analogue of Graver's search algorithm ^(j- I n contrast, 
the purpose of the present paper is to derive more generic 
results for an Hamiltonian-based algorithm perturbed by 
a noise that is described by a stationary gaussian ran- 
dom process. This makes it possible to carry out a fully 
analytical treatment of the tolerance to noise, although 
this is at the price of some (fairly general) assumptions 
on the noise model and of the use of perturbation theory. 

This paper is organized as follows. In Section^ we de- 
scribe our theoretical model of noise based on the Gaus- 
sian Orthogonal Ensemble. In Section IIIII we use per- 
turbation theory to analyze the effect of noise on a time- 
independent Hamiltonian evolution and apply our results 
to the analog analogue of quantum search. In Section llVl 
we consider the tolerance to noise of a time- dependent 
quantum algorithm by adiabatic evolution, and then fo- 



cus on the quantum search by local adiabatic evolution. 
Finally, in Section [V] we conclude by discussing the scal- 
ing of the noise-induced error probability as a function 
of the noise bandwidth. 



II. NOISE MODEL 

Suppose we have an ideal (noiseless) Hamiltonian al- 
gorithm based on the Hamiltonian H (t), 

ihj t \m) = m\m)- (i) 

At the end of the computation (t = T), we obtain the 
state |t/>(T)), which, after measurement, defines the out- 
put of the algorithm. Now, suppose that a perturbation 
eh(t) adds to the ideal Hamiltonian: 

H(t) = H(t) +eh(t). (2) 

Instead of |?/>(T)), we will get at the end of the com- 
putation a different state \ip(T)). The problem in the 
following will be to evaluate the error probability 

Pc rr = l-|WTMT))| 2 (3) 

induced by the perturbation. 

In order to derive analytical results, we will have to 
make some assumptions on the noise term eh(t). First, 
we limit ourselves to a noise of small amplitude e <C 1, 
so that the use of perturbation theory is justified. Sec- 
ond, we assume that, in any basis \<fk) (k =, . . . , N — 1) 
of the ./V-dimensional Hilbert space where the computa- 
tion takes place, the matrix elements of h(t) are normal 
random variables: 

hki(t) = ((fik\h(t)\ipi) 6JV(0,4). (4) 

More specifically, we assume that the matrix h(t) is 
drawn from a Gaussian Orthogonal Ensemble (GOE), 
so that the variance <rjy of its elements is defined by 
°ki = (l + ^fci) ' 2 ! where a 2 is an overall variance (see [Tl| 
for more details on random matrix ensembles). Moreover, 
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any two distinct elements of a GOE matrix are taken as 
independent random variables, 



(h k i(t)h k u>(t)) = 



(M)^(*'.0 and {k,l)^(l',k'). 



(5) 



Even though the above assumptions are not based on a 
specific physical source of noise, they may be justified by 
considering that the noise is generally caused by many in- 
dependent sources of error which, combined together, fi- 
nally result in a random Hamiltonian drawn from a GOE 
as a consequence of the central-limit theorem. 

Furthermore, we assume that the random matrix ele- 
ments hki(t) evolve in time as some stationary random 
process with an autocorrelation function [12j 



R(t) = (h k i(t + T)h k i(t)). 



(6) 



For instance, a very typical noise model we can use is 
a white noise with a high-frequency cut-off lo (see for 
instance 01), which yields 



2 sin o; T 
R ( T ) = a ki ■ 

LOOT 



(7) 



However, to be slightly more general, we will only assume 
later on that the autocorrelation function is of the type 



R(r) = aljM, 



where f{x) verifies f(—x) — f{x), f(x) < /(0) = 1, as 
well as some other regularity conditions (see next sec- 
tions). Thus, we only need to assume that R(t) is a 
function of loqt. 

Finally, as we will be interested in the scaling of the 
perturbed Hamiltonian-based algorithm as a function of 
the size of the problem, N, we need to properly define 
the dependence of the noise term in N. For the scaling 
analysis to be sensible, we must keep a constant signal- 
to-noise ratio as iV increases, that is, the eigenvalues of 
h(t) should scale similarly to those of H(t). As a result 
of Wigner's semi-circular law, we know that the density 
of eigenvalues of GOE matrices for N 1 is given by 



P(E) 



N- 



1 

4a' A tt 







E 2 if \E\ < V4a 2 N 



otherwise. 



(9) 



Therefore, to keep a constant signal-to-noise ratio, we 
have to impose that a 2 = E 2 /AN, where E is of the 
order of the eigenvalues of H(t). 



III. TIME-INDEPENDENT HAMILTONIAN 
EVOLUTION WITH NOISE 

A. Perturbation theory 

Let us study first the simplest case of a time- 
independent Hamiltonian evolution. The solution of the 



ideal Schrodinger equation is 



|^(f)>=£>e- 



\fk), 



(10) 



where \ipt~) and E k are the eigenstates and eigenvalues of 
the ideal Hamiltonian and the amplitudes b k follow from 
the initial conditions. By use of perturbation theory, we 
can study the effect of a small time-dependent perturba- 
tion eh(t) on the ideal Hamiltonian H. Expanding the 
solution of the perturbed equation in the basis formed by 
the solutions of the non-perturbed equation, that is, 



\m) = X>(*k 



\fk) 



(ii) 



and introducing this expression into the Schrodinger 
equation, we get 



b k = 



E 6 



iuikit 



hki(t), 



(12) 



where Lo k i = (E k — E{)/h. Using the same initial state 
as for the ideal evolution, i.e., 6^(0) = b k , we obtain the 
system of equations 



b k (t) 



V f hihy^hu^dh 

, Jo 



(13) 



Using standard perturbation theory (see e.g. |3), this 
may be solved iteratively, building step by step the ex- 
pansion of b k (t) in increasing orders in e. From this so- 
lution, one can derive (an expansion of) the error prob- 
ability p eTI introduced by the perturbation eh(t). As the 
matrix elements of h(t) are random variables, so will p eTI 
be, and we will only have access to its statistics. In par- 
ticular, we will focus on its mean (p eiI ). Using our as- 
sumption that h{t) is a random matrix drawn from a 
GOE, we can show that 



(Perr) = 



= z~ 



k,l 



M 2 (i-M 2 K 



kl 



k^l 



u 2 i 



+0(e 3 ), 



(14) 



where we have introduced the integrals 



"ki 

h 2 



dhdhe^^f^oih-h)), (15) 



which correspond to the coupling between the states \<fk) 
and \ipi) that is effected by the perturbation. Our prob- 
lem now is to evaluate these integrals. We see that they 
only depend on the the noise model via the autocorre- 
lation function f(x) and the high-frequency cut-off loq, 
while they depend on each particular instance of the 
problem via the spectrum of the ideal Hamiltonian or 
the frequencies oj k i (as well as the computation time T). 
Therefore, I kl vary for different instances of a problem. 
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However, we can derive some general expressions, which 
remain valid for a fairly large class of problems. 

First of all, since are integrals over a domain of size 
T 2 and as the amplitude of their integrand is bounded 
by 1, we immediately see that, whatever the values of LOki 
and ujq, we have the upper bound 



< 



'Id 



T 2 



h 2 



(16) 



Furthermore, we note that the 1^ couplings only appear 
between the eigenstates that are initially populated, and 
may therefore be viewed as the interferences caused by 
the noise between these states. As there is in general 
a small and fixed number of eigenstates \ipk) that are 
populated ^ 0) in the algorithm j Eq. I|14|) implies 
that there will be a fixed number of terms contributing 
to the expression of (p crr ). In contrast, the number of 
1^ terms, corresponding the the coupling of the initially 
populated states to all others, will in general grow with 
the dimension N of the Hilbert space. Therefore, the 
scaling of the average error probability (p e rr) will mostly 
depend on the integrals 1^, which is why we now focus 
on these in what follows. By changing the integration 
variables to u = t\ — £2 and v = t\ + we get 



I kl =2 fit J dv J ducos(uj k iu)f(uj u), (17) 

which is the integral of a modulated oscillation. 
For a white noise Q, we get by direct integration 



'Id 



W H 2 LU 



1 - cos(uj k i - uj )T 



OJkl - ^0 
1 - cos(uj k i + LU )T 



+ T Si{uJki - w )T 



(18) 



where Si(ir) is the sine integral function. Depending on 
the value of u>q, we may consider two limiting regimes: 
for a high cut-off frequency uq 3> u>ki , we get 



L kl 



h 2 L 



o 



1 



WW 
UJ 



U! T 



(19) 



while for a low cut-off frequency uj -C u^i , we have 



n 2 



wr. 



01 + ^ 



Wp 
UJkl 



(20) 



Although Eqs. (fTT))) and (|2"0")l are only valid, strictly 
speaking, for a white noise, we obtain similar results for 
a general function f(x). In the high-wo regime, since the 
autocorrelation function R(t) usually tends to zero as 
r increases [i.e., hki(t + r) becomes less and less corre- 
lated with hki (t) for increasing r] , Eq. H19fl follows from 
the approximation cos(ojkiu) f (ujq u ) = 1 + O(u>ki/too)- In 
the low-wo regime, we must integrate a rapidly oscillat- 
ing function over many periods, which is treated in Ap- 
pendix ^ Under very general regularity conditions on 



f(x) 0|, we may use the lemma |2 twice, and finally re- 
cover Eq. H2U|) . It is interesting to note that in the low 
frequency regime, the coupling integral does not depend 
on the computation time T. We will see that, at least for 
the algorithms considered here, this causes a very differ- 
ent behavior of the scaling for this regime, as compared 
to the high-wo regime. 



B. Analog quantum search 

Let us recall the principle of Farhi and Gutmann's ana- 
log quantum search |lOj . Suppose that we may apply an 
oracle Hamiltonian 



H f = E(I-\m)(m\). 



(21) 



to the system, with E denoting an energy scale of the 
system. The problem is to prepare the system in the 
(unknown) solution state |m). In |lOj . it is shown that 
this may be achieved by preparing the system in the uni- 
form superposition of all states 



JV-l 

!*>> = ^ E 1*) 

v fc=0 



and applying the constant Hamiltonian H = Hq 
where 



during a time 



H = E(I-\rM(rl>o\), 



„ nh 1 — 
T = ^VN. 
2E 



(22) 



H f 



(23) 



(24) 



This results in a quadratic speed-up with respect to a 
classical search in an unstructured database of size N. 

In order to study the robustness of this quantum al- 
gorithm against a stationary gaussian noise as defined 
in Sec. |Hj let us first consider the spectrum of the ideal 
Hamiltonian H = Hq + Hf (see Fig. QJ. We assume, for 
simplicity and without loss of generality, that the prob- 
lem admits the solution m = 0. The two lowest eigen- 
values of H, that is Eq = (1 — x)E and E\ = (1 + x)E 
with x = 1/ \/N, are non-degenerate and correspond to 
the ground and first-excited states, 



Wo) 



\<Pi) 



1 + x 



l-x 



N-l 



10} 



10) 



VWTxjf, 



N-l 



E l fc >' ( 26 ) 



2E 



whereas the N — 2 times degenerate eigenvalue E2 
corresponds to the eigenstates 

|^) = -L(|fc) - |1)) k = 2,---N. (27) 
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FIG. 1: Spectrum of the ideal Hamiltonian H — Ho + Hf. 



Expressing \ip(t = 0)) = j^o) i n terms of the eigen- 
statcs \ifk) of the ideal Hamiltonian H, we get 

WO)) = y^ko) - \[^\<Px) (28) 

As a consequence, the instantaneous state of the ideal 
algorithm \ip(t)) is given by Eq. l|10f) with 

- / 1 + x - / 1 — x .„„. 

^ = V— ' 5i = -V^' (29) 

and bk = for k > 2. Only two states are populated dur- 
ing the ideal algorithm, and the average error probability 
(I14|l becomes 

(p err ) = e 2 {(N - 2) [l&ol 2 ^ + Ihfl^] 

+ l&o| 2 |6i| 2 (/oo+^i) + (l&o| 4 + |6i| 4 )/oi 
-2Re [(b* b 1 ) 2 I+] } • (30) 

where we have used the normalization condition |&o| 2 + 
| &i| 2 = 1 and the fact that — Ij~ 2 for I > 2. For this 
algorithm, the bound l(TB|l gives |-fjy(t)| < tt 2 /8, which 
is independent of N. Therefore, only the first term of 
Eq. 1(50) 1 , which represents the coupling of the ground and 
first excited states (the only initially populated states) to 
the N — 2 others, can grow with N and must be taken 
into account in the scaling analysis. Let us focus on this 
term in the two limiting regimes considered above. 

For a noise with a high cut-off frequency, Eq. {TQJ 
yields 

'«=i&°(^( i+ £0 " i = o ' i) (3i) 

which is valid if Hluq 3> E. Clearly, I7„ should be of 
order 1/N for (p e rr) not to grow with N, which imposes 
the condition 

Hluq > Ey/N. (32) 
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FIG. 2: Average error probability (p e rr) (to the second order) 
due to a noise modeled as in Sec. IH] with an autocorrelation 
function f(x) — sin 2/2; for the analog search among iV 
I UN elements. Note that (p CII ) stays very small as long as 
Tiujq/E < 1, at which point it shows a sudden increase. For 
larger values of Kujo, it tends progressively back to a low value. 
The peak at huio ~ E scales as y~N . 



Thus, in this regime, the cut-off frequency of the noise 
must increase as the square root of the size of the prob- 
lem in order to keep a probability of error of constant 
order. In the case of a noise with a low cut-off frequency, 
Eq. JDJ yields 

I k2 = jt°( 1 +ir) ( k = ^ ( 33 ) 

so we see that (p e rr) will not grow with N as long as 
hujQ <C E. Interestingly, this upper bound on the cut-off 
frequency does not depend on the size of the problem. 

We conclude that the influence of noise on the analog 
quantum search algorithm is negligible if the noise varies 
either very slowly or very rapidly with respect to the nat- 
ural time scale of the problem h/E. For the typical case 
of a white noise with high frequency cut-off, the exact in- 
tegration i|18|) shows that for intermediate values of the 
cut-off frequency ujq <~ E/h, the average error probabil- 
ity scales as e 2 ^/N for a given signal-to-noise ratio e (see 
Fig. 0) . This means that there exists a forbidden band 
for the cut-off frequency c^o if we want to keep our algo- 
rithm robust with respect to a constant noise (i.e., a noise 
that scales with N so to keep a constant signal-to-noise 
ratio) . Alternatively, we see that the signal-to- noise ratio 
e should scale as TV -1 / 4 in order to keep the error prob- 
ability (p e rr) constant for a cut-off frequency luq of order 
E/h. In other words, for an increasing problem size N, 
the noise variance cr 2 must decrease faster than E 2 /AN 
for the algorithm to remain immune to noise. Finally, 
let us mention that e ~ TV -1 / 4 coincides with the result 
of Shenvi et al. Q, even though their noise model was 
less general than ours since they only considered an error 
in the magnitude of the oracle Hamiltonian, modeled as 
a Markovian stochastic variable with Gaussian distribu- 
tion. 
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IV. ADIABATIC EVOLUTION WITH NOISE 

A. Adiabatic approximation 

Let us recall the adiabatic approximation, which is at 
the basis of the quantum algorithms by adiabatic evolu- 
tion. Qualitatively speaking, the idea is as follows: if a 
quantum system is prepared in its ground state and its 
Hamiltonian varies "slowly enough" , it remains in a state 
close to the instantaneous ground state of the Hamilto- 
nian at any time. To be more precise, let us consider the 
Schrodinger equation for a time-dependent Hamiltonian 
(see Q for details), 



%h-w))=H(t)m)). 



(34) 



To solve this equation, we express its solution |$(t)) m 
the basis formed by the instantaneous eigenstates \ipk(t)) 
of the Hamiltonian Hit), 



\m) 



E< 

k 



hi^e-^o^ 11 ^^)), (35) 



where E^it) are the corresponding instantaneous eigen- 
values of Hit). By inserting this expression into Eq. (|34p. 
we find the system of differential equations 

m = ^Htyti*^ ^{f iff • (36) 

If the quantum system is initially in its ground state 
|V>(G)) = |yo(0))i these equations can be integrated, giv- 
ing 



h(t) = 5 ( 



(I/,- 



(37) 

J^Jo Ek(ti) - hi(ti) 

As in perturbation theory, these equations may be solved 
iteratively, which gives after one iteration 



4"./.= / - iI ^^ {V E^^ dtl (38) 



for k 7^ 0. Now if the variation dH/dt of the Hamiltonian 
is slow enough, or more specifically if 



\Mt)\ 



. |(y fc (t)|#|y>o(f))| 
' (E k (t) - E Q it)f 



< 6 k < 1, (39) 



and under suitable regularity conditions, we may inte- 
grate Eq. I)38|l by parts as shown in Appendix EI which 
yields 



b£\t) = -i Wit^jti 1 ^ dt '} + 0(5l) 



(40) 



We see that, at the first order, the amplitudes |&fe(t)| are 
bounded by 25^- This first-order (so-called adiabatic) 
approximation is acceptable if |&o(i)| 2 stays close to 1, 
that is, if the probability pit) = Ylk^o l^fcWI 2 of hopping 
to any excited state remains small. In other words, if the 
adiabatic condition 



45>u P L4 fe (t')| 2 <S 2 , 



(41) 



is satisfied, then pit) < 5 2 <C 1, with J « 1 being a 
"slowness" parameter. 

Now, suppose that a time-dependent perturbation 
eh(t) adds to the ideal Hamiltonian Hit). We again ex- 
press the solution of the perturbed Schrodinger equation 
in the basis formed by the instantaneous eigenstates of 
Hit), 

\m) = xx*)e- i; ° ^ dti \Mt)) (42) 



which transforms Eq. (J23) into 

bkit)=S ok (43) 

.f«wf«i^ <^(*0lglw(*0> 



l^k 
.£ 



E / 



E k ih) - Efih) 



-dti 



4E / ^C*!)^^ 1 ^*' 1 ^^)^^)!^^))*!- 



We may again solve this system of equations iteratively, 
which gives after one iteration 



(i) 



(44) 



for k ^ 0. As before, this first-order approxima- 
tion remains valid provided that the probability pit) = 
Sfc^o °f hopping to any excited state remains 

small. 

Let us now evaluate the average error probability at 
the end of the evolution t = T using the same model 
as before for the perturbation hit). Defining the error 
probability as p crr = piT), we have 



e 



i/* 1 u> k0 (t')dt' 



(Perr) = Pen + ^2 E / / dt l dt ^ J ' 2 

X {{Mti)\Hti)\Mti)){<pkit2)\hit 2 )\<poit2))) 

+ Oie 3 ), (45) 

where p err — piT) is the error probability of the ideal 
adiabatic evolution. Let us note that \ipkit\)) ^ ^kit-ij) 
in general, so that we do not immediately recover the 
autocorrelation function of one particular matrix element 
of hit). However, as the different matrix elements of hit) 
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are independent in a particular basis, we have 
= ((<Pk(h)\<Pk(ti)) (<po(h)\<po(h)) 

+ (<Pk(t2)\<Po(tl))(<Pk(tl)\<Po(t2))) 

x*l f(u, (h-t 2 )). (46) 

If \ipk(t)) varies sufficiently smoothly for < t < T, 
then the first factor is of order 1 - 0{{h ~ t 2 ) 2 /T 2 ). 
Thus, we may approximate it by 1 as long as ujqT 3> 
1, that is, if the noise varies quickly compared to the 
adiabatic evolution. In that case, we get for the average 
error probability 

(pcrr) = Perr + S 2 £ I w + 0((6 + £ ) 3 ), (47) 
fc^O 

where the integrals 

= JJ^dhdhe 1 ^^ ^ dt 'f(LJ (h-t 2 )), (48) 

represent the coupling of the ground state to the excited 
states induced by the perturbation. The adiabatic con- 
dition generalizes in the case of such a perturbation to 




4sup|,4 fe (t)| 5 



s 2 I, 



k() 



<S 2 



(49) 



Similarly to the case of the perturbed time-independent 
Hamiltonian evolution, the effect of the perturbation on 
the adiabatic evolution mainly depends on the coupling 
integrals I^ Q1 which are bounded by 



I4" I< 



doll 

h 2 



(50) 



and can be approximated as 



'fcO 



l k0 



H 2 G 



o 



uj T 



LUq > LU 



A-0 



^fcO — 



'fcO 



rO 1 



^0 
,min 
J k0 



u < u^a 



(51) 



(52) 



in the limiting regimes of high or low cut-off frequency loq, 
respectively, where u>$ a < u; k0 (t) < u™ ax for t G [0, T], 



B. Adiabatic quantum search 

The principle of the adiabatic quantum search is 
to apply the Hamiltonian Hq = E(I — |V'o)('0o|) to a 
system prepared in its ground state | •i/ , o) an d then to 
progressively switch the Hamiltonian Hq to the Hamil- 
tonian Hf = E(I — |m)(m|), where m is the solution of 



the search problem. If this switch is done slowly enough, 
the system will stay in the instantaneous ground state 
of the Hamiltonian and thus end up in the ground state 
of Hf, i.e., the solution state |m). The instantaneous 
Hamiltonian is chosen as 



H(s) = (l-s)H + sH f , 



(53) 



where s = s(t) is an evolution function which must be 
optimized so as to reduce the computation time while 
respecting the adiabatic condition l|4*T)> (see 0, 0] for 
details). In this case, this condition may be rewritten as 



ds_ < S (E 1 {t) - E (t)) 2 
dt ~ 2h\(ip 1 {t)\Hf-H \ip. (t))\' 



(54) 



Without loss of generality, we may once again suppose 
that Tji — 0. The instantaneous eigenstates of H (s) are 



where 



Ms)) = 

M»)> = 

\<Pk(s)) = 

E (s) = 

E 1 (s) = 

E k (s) = 



NjE^s) - s)\^ ) + s\0) 
V£i(s) 2 + (iV-l)(i?i( S )-s) 2 

VN(E (s) - s)\^ ) + s\0) 
^/Eo( S ) 2 + (N-l)(E Q ( S )- S ) 2 



E 
2~ 

E 
~2 

E 



N-l , 
I — \ I — 4 — s(l — s) 



N-l 

1 - 4^rr-*(l - «) 



N 



k>2 



(55) 

(56) 
(57) 

(58) 

(59) 
(60) 



are the instantaneous eigenvalues of H(s) (see Fig- EJl - 
Since \ \Hf — Hq\\ < E, taking an evolution function s(t) 
that satisfies 



ds 
dt 



2hE 



(E^-EoW 



(61) 



complies with the adiabatic condition l|54|l . and then 
leads to a computation time 



5 E 



(62) 



Consider that some noise, modeled again as a sta- 
tionary gaussian random process, perturbs the evolution. 
Equation (|47JI then reads 

(Pexr) < Pcrr + S 2 [/f + (N - 2)1'] + 0((S + ef), (63) 

where < 7r 2 /(64<5 2 ) as a result of Eq. (JSUJl. Let us 
emphasize that, while it was only the excitation to the 
first excited state that was critical for the ideal adiabatic 
algorithm, in this case it is the coupling of the ground 
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FIG. 3: Instantaneous eigenvalues of H(s) for N — 32. 



state to all excited states that could make the algorithm 
fail, since their number grows as the size N of the prob- 
lem. Moreover, this bound already suggests that the cou- 
pling integrals Ir Q - and therefore the error probability 
- could increase when the evolution slows down (<5 de- 
creases) which means there must be a compromise be- 
tween a slow evolution, very close to perfect adiabaticity, 
and a fast evolution, more robust to noise. 

As before, let us consider the two limiting regimes of a 
high or a low cut-off frequency uiq. In the case of a high 
cut-off frequency (Huj ^> E), Eq. (|ST)) yields 

J ko = 7^0 (-L= (l + ^L\) V k f 0, (64) 

exactly as for the analog quantum search except for the 
factor 1/5. The latter factor shows that in order to keep 
the algorithm robust to noise, the cut-off frequency has to 
increase not only as the the size of the database TV grows 
(just as for the analog quantum search), but also as the 
evolution slows down (5 decreases). More precisely, we 
see that the perturbed adiabatic condition l|49|l is satis- 
fied only if 

E 2 - 

hu) > -~eVn (65) 

(compare with Eq. When the cut-off frequency 

becomes very low (Hcoq <C E), Eq. (|52|l implies that the 
coupling integrals behave as 

for all excited states except the first one (i.e., for k > 2). 
For the first excited state (k — 1), we have fku^ n ~ 
E/y/N so Eq. does not yield a useful result. In- 
stead, we simply use the general bound < n 2 / (646 2 ) . 
Therefore, the adiabatic condition is satisfied here as long 
as Hu)q <C E (just as for the analog quantum search), but 
also if e <C S. 



In summary, we recover essentially the same effects for 
the adiabatic quantum search as for the analog quantum 
search, that is, the influence of noise becomes negligi- 
ble only in the case of a very high or a very low cut-off 
frequency u>q, apart from the influence of the slowness 
parameter 8. Regarding this latter parameter, we see 
that while decreasing <5 gets the ideal evolution closer to 
adiabaticity and therefore reduces the error probability 
without noise, in the presence of noise it imposes that 
e decreases - i.e., that the signal-to-noise ratio increases 
- in both regimes of a high or low cut-off frequency (or 
that the high cut-off frequency increases as l/c> 3 ). 



V. CONCLUSION 

We have studied the resistance of Hamiltonian quan- 
tum algorithms (including adiabatic algorithms) to a 
noise that is modeled as a random matrix whose elements 
are stationary gaussian random processes within a fixed 
bandwidth. This statistical noise model is generic, and 
should therefore make our analysis valid over a large class 
of physical systems, regardless of the exact origin of the 
added noise. Another main advantage of this noise model 
is that it makes it possible to perform a fully analytical 
scaling analysis. Our general result is that the Hamil- 
tonian algorithms are resistant to noise (i.e., the error 
probability (p e rr) does not increase with increasing prob- 
lem sizes N) as long as the cut-off frequency of the noise 
is either very high or very low with respect to the inverse 
of the characteristic time-scale of the system E/H. Asides 
from the influence of the slowness parameter S in the case 
of the adiabatic algorithms, this resistance is essentially 
similar for adiabatic and time-independent Hamiltonian 
algorithms. Our results are in good agreement with the 
numerical study of Childs et al. inlg . They even corrob- 
orate the results of Shenvi et al. 9] , although their noise 
model was rather different, which supports the idea that 
using random matrix theory provides a rather general 
description of noise. 

Roughly speaking, the two limiting regimes of high or 
low cut-off frequencies can be understood in the follow- 
ing way. If the frequency components of the noise are 
much below the inverse of the characteristic time-scale 
of the system, it is intuitively clear that the noise can- 
not effect transitions to undesired states. On the con- 
trary, if the noise spectrum spreads over a band which is 
much broader than the inverse of the characteristic time- 
scale of the system, then the noise spectral density is low 
around the frequencies that effect undesired transitions. 
In the intermediate region, we found that the error prob- 
ability unfortunately scales as \/~N [2£|, which implies 
that some error correction is needed to make Hamilto- 
nian algorithms scalable. This last point is particularly 
important as it is plausible that the source of noise oc- 
curring in a physical system typically varies on a time 
scale comparable to the natural time scale of the system, 
so that the less favorable regime (fouiQ ~ E) may be the 
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most common situation. 

Coming back to the two limiting regimes, let us notice 
that the high cut-off frequency increases towards higher 
frequencies when TV raises [see Eqs. Ij32|l and Ij65|l ]. Con- 
sequently, the Hamiltonian algorithms will in practice not 
be scalable in the high cut-off frequency regime too, since 
the noise should spread over an arbitrarily large spectrum 
to keep the error probability low. In this case, some kind 
of error correction should also be implemented if the size 
of the problem becomes too large. The case of a low cut- 
off frequency, however, is more favorable. Indeed, the sit- 
uation is quite different here as the error probability stays 
small as soon as the cut-off frequency is lower than some 
fixed value, even when the size of the problem increases. 
This fault tolerance may be explained by the fact that 
the spectral density of noise does not contain frequencies 
close to resonances, and thus will not efficiently couple 
different eigenstates of the ideal Hamiltonian. 

It should be emphasized that it is not the possible 
excitation to one particular state that makes the algo- 
rithm fail, but the fact the dimension of the Hilbcrt 
space increases with the problem size, and hence the 
number of states that could be accidentally populated 
as well. This means, in the case of adiabatic computa- 
tion, that even if the gap between the ground and first 
excited states decreases, the algorithm may remain ro- 
bust to a noise with a low cut-off frequency (even if this 
frequency remains constant) as long as the gap between 
the ground states and the other excited states remain 
lower bounded. Therefore, the algorithm would remain 
scalable in the case of a low cut-off frequency as long 
as the natural frequencies of the ideal Hamiltonian are 
much larger than the frequencies contained in the noise. 
Of course, throughout this analysis, we always made the 
assumption that the signal-to-noise ratio remains essen- 
tially constant when the size of the Hilbert space where 
the computation takes place becomes large, which may 
practically not be the case. Thus, even in this low cut- 
off frequency regime, it may be necessary to devise error 
correction techniques for Hamiltonian - and in particular 
adiabatic - quantum algorithms. 

Note: A few days ago, a related paper appeared on the 
quant-ph preprint server, which also considers the use of 
random matrix theory in adiabatic quantum computing 
but with a distinct goal, namely to analyze the spectral 
statistics of the Hamiltonian over a large class of prob- 
lems along an adiabatic (but noiseless) evolution 17]. 
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APPENDIX A 

In this appendix, we give a useful tool to evaluate in- 
tegrals of an oscillating function such as 



dxF{x)e lux . (Al) 

The basic idea relies on Riemann-Lebesgue's lemma: 

Lemma 1 (Riemann-Lebesgue) // F(x) is an inte- 
grable function on [a,b], then 



lim 

U) — >oc 



dxF(x)e iuJX = 0. 



This lemma suggests that the integral (|Alll will be rel- 
atively small if lu is sufficiently large. The purpose of this 
appendix is to quantify this idea. 

First of all, as long as F(x) is differentiable on [a, b], 
we may integrate (IA1|) by parts: 



dxF(x)e i 



[F(x)e 



dF 

dx—{x)e lux 1 
dx 



where [/(x)] a = /(b) — /(a), and, using this last equation 
iteratively, we show that for an Af-timcs differentiable 
function F(x) on [a, b], 



N-l 



dxF(x)e i ' 



+ 







'd n F 


(= 




dx n 


JV 


r b d N F 


) , 


f dx- 


dx N ( 




a 



{x)e lux . (A2) 



The order of the error introduced by neglecting the last 
term may be evaluated as follows: 



N 



dx . N [x)e 



dx 



d N F 



dx 



N 



(X) 



We see that the accuracy of this approximation in- 
creases with the oscillation frequency to. Moreover, if 
(1 / lo n )d N F / dx N — ► for N — ► oo, this error approaches 
zero as N increases and we prove the following lemma: 

Lemma 2 Let the function F[x) be infinitely differen- 
tiable on [a,b]. If 



1 d N F 



dx N 



to 



N 



(X) 



for some real u>, then 
dxF{x)e iux = - 



n=0 



N- 



V x 



n+l 



d n F 

dx n 



While this result is helpful to study a time- independent 
Hamiltonian evolution, in the case of an adiabatic evo- 
lution the typical frequencies become time-dependent. 
However, using the same method, we easily generalize 
this lemma to the case of a varying frequency lo(x). 
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Lemma 3 Let the function F{x) be infinitely differen- 
tiable on [a, b]. If 



1 d N F 



(,•) 



V .x 



Lu(x) N dx N y ' JV^oo 

for some real differentiable function lo(x) on [a,b], then 



dxF{x)e l ^ u]{x ' )dx ' 



£ 

71 = 



n+1 



d n F 



u)(x) J dx n 



(x)e i f° u ^ x '^ dx ' 



a ^ dx W 35 ) 



n+1 



d n F 

dx n 



(x)e l f« ^ x "^ dx ' 
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